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Abstract. We extend the Langevin approach to a class of driving noises 
whose generating processes have independent increments with super-heavy-tailed 
distributions. The time-dependent generalized Fokker-Planck equation that corres- 
ponds to the first-order Langevin equation driven by such a noise is derived and solved 
& ' exactly. This noise generates two probabilistic states of the system, survived and 

absorbed, that are equivalent to those for a classical particle in an absorbing medium. 
' The connection between the rate of absorption and the super-heavy-tailed distribution 

of the increments is established analytically. A numerical scheme for the simulation of 
the Langevin equation with super-heavy-tailed noise is developed and used to verify 
i our theoretical results. 
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1. Introduction 

At the beginning of the last century, Paul Langevin, a prominent French physicist, 
proposed an alternative description of Brownian motion [TJ. In contrast to earlier 
Einstein's approach [2j, which deals with the Fokker-Planck equation for the proba- 
bility density of an ensemble of independent Brownian particles, Langevin used the 
equation of motion for a single Brownian particle, the so-called Langevin equation. In 
fact, applying Newton's second law with a random force arising from the surrounding 
medium, he introduced a first example of stochastic differential equations. Now the 
Langevin approach is one of the most effective and widely used tools for studying the 
effects of a fluctuating environment [3]. 

In physics, many other stochastic equations that account for the influence of 
the environment are also often termed "Langevin equations". Among them, the 
dimensionless first-order (overdamped) Langevin equation 

x(t) = f(x(t),t)+m (i) 

is the simplest and at the same time most important. It describes noise effects in very 
different systems, but to be specific we will call x{t) a particle position [x(0) = 0], f(x, t) 
a deterministic force field, and £(t) a random force (noise) resulting from a fluctuating 
environment. The Langevin approach is particularly useful when the noise can be 
approximated by the time derivative of the noise generating process rj{t), i.e., random 
process with stationary and independent increments Ai](t) = r/(t + r) —i](t). The reason 
is that in this case the solution of equation ([1]) is a Markov process whose future behavior 
is determined only by the present state, see e.g. [HE]. 

The main statistical characteristic of the particle position x(t), its probability 
density P(x,t), depends fundamentally on the probability density p(Ar], r) of the 
increments Ai](t) for a sufficiently small time increment r. On a phenomeno logical 
level, p(Ar], t) must be chosen to obtain the best representation of the environmental 
fluctuations. If these fluctuations are the result of many independent random variables 
with zero mean and finite variance, then, according to the central limit theorem [7J, 
p(Ar], r) can be approximated by a Gaussian probability density with zero mean and 
a variance proportional to r. Accordingly, r](t) becomes a Wiener process and 
reduces to Gaussian white noise. As is well known, in this case the process x(t) 
is continuous and P(x, t) satisfies the ordinary Fokker-Planck equation [3] [6] which 
describes a huge variety of noise phenomena. If the above mentioned random variables 
have no finite variance then, as the generalized central limit theorem suggests [7J, 
p(Ai], t) can be chosen in the form of a Levy stable probability density. This choice of 
p(An, t) corresponds to Levy stable noise and leads to the fractional Fokker-Planck 
equation for P(x, t) [8] [12]. The Langevin equation ([1]) driven by Levy stable noise and 
the corresponding fractional Fokker-Planck equation describe the so-called Levy flights, 
i.e., random processes exhibiting rare but large jumps [T3] [T6]. 

Recently, it has been shown [17] that the ordinary and fractional Fokker-Planck 
equations hold also for the probability densities represented in the form p(Ai], r) = 
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p(Ar] I a(r)) / a(r) , where a(r) is a scale function that tends to zero as r — > 0, and 
p(y) is a probability density satisfying the condition lim^o p(y I e ) I e — 8(y) with 5{y) 
being the Dirac delta function. Specifically, the ordinary equation holds for a class of 
probability densities p(y) with finite variance, and the fractional one for a class of heavy- 
tailed p(y), i.e., probability densities with power-law tails and infinite second moment. 
Although these classes cover the most important cases, they do not exhaust all possible 
probability densities: There also exists a class of super-heavy-tailed probability densities 
p(y) which have never been used in the Langevin approach. According to the definition, 
these densities decay so slowly that all fractional moments M 7 = Jf^ dy |y| 7 p(y) with 
7 > are infinite (at 7 = the normalization of p(y) implies Mo = 1). Because 
of this "unphysical" property, the super-heavy-tailed distributions, as well as the 
heavy-tailed ones, are not directly related to "physical" random variables having finite 
variance. Nevertheless, they are often used to model some specific properties of physical 
systems. Moreover, the heavy-tailed distributions have become a common tool for 
studying anomalous physical phenomena, see e.g. p2]-[T5] and [18]-|20j. Although 
the super-heavy-tailed distributions have not yet received much attention, the examples 
of extremely slow diffusion [211 E2] demonstrate the usefulness of these distributions for 
describing systems with highly anomalous behavior. 

The purpose of this paper is twofold. First, we wish to extend the Langevin equation 
method to a class of super-heavy-tailed noises, i.e. noises arising from the super-heavy- 
tailed distributions. The importance of this problem is that these noises, together 
with the previously known, exhaust all possible noises associated with the probability 
density p(y) and, as a consequence, enhance the capability of this method. Second, 
we intend to show that super-heavy-tailed noises generate two probabilistic states of 
a particle, survived and absorbed. The existence of these states makes it possible 
to model the particle dynamics, which is randomly interrupted by the transition of 
a particle to a qualitatively new state, within the Langevin approach. This transition 
can occur in a number of physical processes, including the processes of annihilation, 
evaporation and absorption. In this paper we demonstrate the applicability of the 
Langevin equation ([1]) driven by super-heavy-tailed noise for describing the particle 
dynamics, both deterministic and random, in an absorbing medium. 

2. General results 

We start our analysis with the generalized Fokker-Planck equation [23] 

^P(x, t) = ~f(x, t)P(x, t) + T- x {P k (t) In S k } (2) 

which corresponds to the Langevin equation (JTJ. Here F{u(x)} = u k = /f^, dxe~ lkx u{x) 
and J r_1 {nfc} = u(x) = (1/2-71") Jf^ dk e lkx u k denote the direct and inverse Fourier 
transforms, respectively, P k {t) = J^° 00 dxe~ 1,kx P(x,t) is the characteristic function of 
the particle position x(t), and S k is the characteristic function of the noise generating 
process r)(t) at t — 1. It is assumed that the solution of equation (j2j) satisfies the initial, 
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P(x, 0) = S(x), and normalization, J^ OQ dxP(x,t) = 1, conditions. The characteristic 
function S k depends on the Fourier transform of the transition probability density 
p(Ai],t) over the variable Arj. In particular, if p(Ar), r) = p{Ar] / a{r)) / a{r) then 

mSfc = lim-[p fco ( T ) - 1] (3) 

with Pfe a ( T ) = J^dy e~ lka( - T ^ y p(y). Equation reproduces all known forms of the 
Fokker-Planck equation associated with the Langevin equation (CQ) and is valid for all 
probability densities p(y) [21]. It should be noted that because the scale function a(r) 
tends to zero as r — > 0, the dependence of In Sfc on k is determined by the asymptotic 
behavior of p{y) at y — > ±oo. 

Next we consider a class of symmetric super-heavy-tailed probability densities p(y) 
whose asymptotic behavior at y — > oo is given by 

p{y) ~ -%), (4) 

y 

where h(y) is a slowly varying positive function. The term "slowly varying" means 
that the asymptotic relation h(fiy) ~ h(y), i.e. lim^^ h{\xy )/h(y) = 1, holds for all 
/i > 0. Since the probability density p(y) is normalized, the function h(y) must tend to 
zero more rapidly than l/\ny. On the other hand, since y u h{y) — > oo as y — > oo for 
all v > [25], must decay more slowly than any positive power of 1/y, yielding 
M-y = oo (7 > 0). An illustrative example of the probability density p{y) satisfying 
these conditions is given in equation ffT4"]) below. 

In order to find In for the examined class of probability densities, we first rewrite 
([3]) in the form lnS^ = — 21im T ^ ^(l/l^l a ( r ))/ r with 

Y(X)= I dy[l-cos(y/X)]p(y) 

POO 

dxsinx / dyp(y) (5) 

Using either of these representations and the asymptotic formula (T4]), one can easily 
show that the function Y(X), which according to (J5j) tends to zero as A — > 00, is slowly 
varying, i.e. Y(l/\k\a(r)) ~ Y(l/a(r)) if k 7^ and r — > 0. Finally, taking into account 
that In So = 0, we obtain the general form of In Sk in the case of super-heavy-tailed p(y): 

kiS k = -q(l- 8 k0 ), (6) 

where 5 k o is the Kronecker delta which equals 1 if k — and otherwise, and 

g = lim-y(-M. (7) 

r^or \a(r)J 

This limit depends on both the asymptotic behavior of Y(X) at A — > 00 and the behavior 
of the scale function a(r) at r — > 0. Since the former is controlled by a given function 
h(y), the nontrivial action of super-heavy-tailed noise, which is characterized by the 
condition < q < 00, occurs only at an appropriate choice of a(r). We note that while 
a(r) oc T 1 / 2 for all p(y) with finite variance and a(r) oc r 1//a almost for all heavy-tailed 
p(y), where a is an index of stability associated with the corresponding p(y) [17], there 
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is no universal dependence of a(r) on r for different super-heavy-tailed p(y). Moreover, 
because of the function Y(X) is slowly varying, the scale function a(r) is in general not 
unique even for given p(y) and q. We stress, however, that according to (jEJ) the influence 
of super-heavy-tailed noises is fully accounted by the parameter q. Therefore, this lack 
of uniqueness is of no importance and both the Langevin and Fokker-Planck equations 
remain well defined (see also the next section). 

To avoid possible problems with the interpretation of the Fokker-Planck equation 
(j2J at the condition (jSJ), let us temporarily replace the Kronecker delta by a less singular 
function A k (A*. = 1 if \k\ < k and A k = if \k\ > k) which reduces to S k Q in the limit 
k — y 0. In this case, taking In S k = — q(l — A k ), we can write the solution of equation (J2J) 
in the form P(x, t) = V(x, t) +A(x, t), where the terms V(x, t) and A(x, t) are governed 
by the equations 

^V(x, t) = ~f{x, t)V(x, t) - qV(x, t), (8) 

—A(x, t) = -faf(x, t)A(x, t) - qA{x, t) + qF-^Pkit) A k } (9) 

with the initial conditions V(x, 0) = 5(x) and A(x, 0) = and the normalization 
condition Po(t) = Vo(t) + Ao(t) = 1. Assuming that z(t) is the solution of the 
deterministic equation z(t) = f(z(t),t) satisfying the initial condition z(0) = 0, the 
solution of equation (jSJ) can be written as 

V{x,t) = e- qt 5[x -z(t)}. (10) 

By integrating V(x, t) over all x one finds that this solution is not normalized to 
unity: Vo(t) = e~ qt . Using this result and the relation Vo(t) + Ao(t) = 1 we obtain 
Ao(t) = 1 — e~ qt . It should be noted that the function V(x,t) and the normalization 
condition for A(x, t), Ao(t) = 1 — e~ qt , do not depend on k, i.e. they hold also at k — > 
when A fc = 5 k0 . 

In contrast, it is clear from the equation 

t) = - t)A {x , t) - ^ t) + *y!: (t y «-' 

7T J-oo x — y 

(A(x, 0) = 0), which follows from ([9]), (TTOT) and the relation 

^{Pk^Ak} = 1-1" dke ikx [V k (t)+A k (t)], (12) 

that the function A(x, t) depends on k. It is not difficult to verify that equation ( TTTT) 
provides a correct normalization for A(x, t). Indeed, the integration of both sides of this 
equation over all x and the use of the formula /f^ dx s\ti(kx)/x = n lead to the equation 
dAo(t)/dt = qe" qt whose solution satisfying the initial condition *4o(0) = reproduces 
the above result: Ao(i) = 1 — e~ qt . We note, however, that because of replacing S k0 by 
A k , equation ( fTTj) with a fixed k describes the influence of super- heavy-tailed noises in 
an approximate way. Moreover, this approximation is valid if \x — z(t) \ is much less than 
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the characteristic length scale 1/ k (in the opposite case the behavior of A(x, t) may be 
unphysical with A(x,t) < for some x and t). 

Assuming that \x — z(t)\ <C 1/ k and k is so small that Vk(t) + Akit) m 1, we can 
reduce equation (fTTj) to the form 

—A(x, t) = -fofix, t)A(x, t) - qA(x, t) + ^. (13) 

As is well known (see e.g. [26J ) , its general solution is given by the sum of a particular 
solution of this equation, which is proportional to k, and the general solution of the 
related homogeneous equation (when k = 0). Since A(x, 0) = 0, the latter is also 
proved to be proportional to k. Thus, the function A(x,t) at \x — z(t)\ C 1/k tends 
to zero linearly with k. In particular, if the external force depends only on time, i.e. 
f(x,t) = /(£), then z{t) = Jq dt' '/(£') and the general solution of equation (1T3"|) is given 
by A(x, t) = e~ qt $(a; — z(t)) + k/tt, where $(x) is an arbitrary function. This function 
can be determined from the initial condition A(x, 0) = 0, $(a;) = —k/tt, yielding 
A(x,t) = k(1 - e~ qt )/ir as k — \ 0. 

The above results show that under the action of super-heavy-tailed noise two 
probabilistic states of a particle appear. The first one is described by the probability 
density V(x,t) and is realized with probability Vo(t). Since the noise does not affect 
the particle trajectory z(t), see ( |T0|) . we refer to this state as the surviving state. The 
second state, which is associated with the probability density A(x, t) (we recall that 
A(x, t) — > as k — > 0), is realized with probability ^4o(^)- The transition to this state 
implies that a particle jumps to infinity and it is excluded from future consideration. It 
is therefore reasonable to call this state absorbing. To avoid any confusion, we emphasize 
that the act of absorption is considered here as the transition of a particle into the state 
in which the probability to find this particle in any finite interval equals zero. Thus, 
super-heavy-tailed noise plays the role of an absorbing medium characterized by the rate 
of absorption q. Accordingly, the Langevin equation (CQ) driven by this noise describes 
the overdamped motion of a particle in an absorbing medium. 

As equation ( flOi) shows, before absorption, which is random in time, the motion is 
deterministic and occurs under a force field f(x, t). The random motion of a particle in 
an absorbing medium can also be described within the Langevin approach. To this end, 
we represent £(t) in the Langevin equation (TjQ) as a sum of two independent noises: (i) 
super-heavy-tailed noise which models an absorbing medium, and (ii) noise £2^), 

which induces the random motion of a particle. Calculations similar to those described 
above lead to the result V(x,t) = e~ gt W(x,t), where W(x, t) denotes the normalized to 
unity probability density associated with the random motion of a particle. In particular, 
if £2^) is Gaussian white noise then W (x, t) is the solution of the ordinary Fokker-Planck 
equation, and if £2^) is Levy stable noise then W(x,t) is the solution of the fractional 
Fokker-Planck equation. In these cases, the Langevin equation describes Brownian 
motion and Levy flights, respectively, in an absorbing medium. 
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3. Illustrative example 

In order to verify numerically that super-heavy-tailed noise acts as an absorbing medium, 
we first specify the super-heavy-tailed probability density p(y) by choosing 

P(v) = , 7 | il^/ , I i\ ? ( 14 ) 
2(s + \y\)\n (s+ \y\) 

where s(> 1) is the density parameter. This symmetric probability density is unimodal 
with the peak located at y = and has no characteristic scale. The behavior of p(y) 
at y — > + is characterized by the conditions p(0) = l/(2slns) and dp(y) / dy\ y ^ + = 
— (1 + 0.5 In s) j is In s) 2 , and at y — > oo by the asymptotic formula p(y) ~ In s/(2y In 2 y) 
[this implies that h(y) ~ In s/(2 In 2 y)]. For this probability density the slowly varying 
function fl5]), after calculating the inner integral, takes the form 

Y(X) = — f^-^— . (15) 
K J 2 Jo \n{s + Xx) v ; 

To find the leading term of Y(X) at A — > oo, which determines the rate of absorption 
([7]), we use the identity 

1 r 1 du 1 

+ rTTTT 16) 



lnt> Jo v u vlnv 
with v = s + Xx. Substitution of this identity into f|T5|) yields 

Y(x) = ^[ 1 d U rdx- sinx 



2 Jo Jo (s + Xx) u 

Ins f°° sin a; . . 

+ ~Jo X {s + Xx) ln(s + Xx)' ^ ' 

Then, taking into account that J °° dx sin(x)/x u = T(l — u) sin[7r(l — u)/2], where T(x) 
is the gamma function, 

1 f°° sinx f 1 r(l — u) n(l-u) 
du / dx- : — — ~ / du : sin 



o Jo (s + Xx) u Jo X u 2 



'^e^--^- (A^oo) (18) 
In A 



^ dx 7 x VTT : — r = o( T ) (A ^ oo), (19) 



and 

/o (s + Ax) ln(s + Ax) V A 

we arrive to the desired asymptotic formula Y(X) ~ Ins (21nA) _1 (A — > oo). 

As it follows from the definition of the rate of absorbtion (JTj), in the considered case 
it can be represented in the form 

?=- (20) 
r 



with 

f = ^ , , , , , 

t->o \a(r) 



limlnl . (21) 
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Depending on the behavior of the scale function a(r) at r — > 0, the parameter r can take 
the values from to oo. If a(r) approaches zero so fast that r = oo then super- heavy- 
tailed noise does not affect the particle dynamics at all, i.e. the particle absorption is 
absent (q = 0). In contrast, if a(r) tends to zero in such a way that r = then the noise 
action is so strong that the transition to the absorbing state happens at t — + , i.e. a 
particle is absorbed immediately (q = oo). A non-trivial action of super-heavy-tailed 
noise occurs at < r < oo. According to (|2"T|) . in this case a(r) ~ b(r)e~ r / T , where 
b{r) is a positive function satisfying the condition [6(t)] t — > 1 as r — > 0. An explicit 
form of this function is of no importance (and so we can choose a(r) = e~ r ^ T ) because 
at a(r) ~ b{r)e~ r l T and lim r _ i ,o[&( r )] T = 1 the rate of absorption q does not depend on 
6(r): 

In s In s In s 

lim^o ln[e r / r /6(r)] T r — lim r ^ ln[6(r)] r r 
Thus, assuming for simplicity that f(x,t) = 0, we obtain for the examined case the 
probability density of the particle position in the surviving (non-absorbing) state, 

V(x,t) = s- t/r 6(x), (23) 

and the probability of absorption, Ao(t) = 1 — s" 1 ^ . 

It is worthy to note that while the probability density ffT4l) with super-heavy 
tails depends on the single parameter s, the corresponding super- heavy-tailed noise 
is characterized also by the additional parameter r. The parameter r relates to the scale 
function a(r) and at a fixed s it controls the noise intensity, namely, the smaller the 
parameter r is the larger is this intensity. The dependence of the scale function on the 
parameter which characterizes the noise intensity is not surprising and occurs for all 
probability densities p(y) [21] . If the parameters s and r are chosen then super-heavy- 
tailed noise associated with the probability density (TH|) is completely determined and 
both the Langevin and Fokker-Planck equations are well defined. 



4. Numerical simulations 



We have verified the analytical results obtained in the previous section by numerically 
studying the statistics of the particle positions. Because of the Langevin equation has 
never before been simulated in the case of super-heavy-tailed noises, we briefly describe 
the procedure for finding the particle position x(nr) = YJj=\ ^ 1 l{j T ) at t = nr. The first 
step is calculating the quantity pj = a(r)(s 1 ^ 1-2 ^ — s) with a(r) = e~ r l T and Zj being a 
random number uniformly distributed in the interval [0, 1]. Since the positive quantities 
Pj are distributed with the probability density Ins {a(r)[s + p/a(r)] ln 2 [s + p/a(r)]} _1 , 
the increments Arj^jr), whose probability density p(Ar] / 'a(r)) / a(r) is determined by 
dHJ), can be represented in the form Ar/(jr) = pj/2, where (±)j = +1 or —1 with 
probability 1/2 each. Finally, calculating At)(Jt) n times, we obtain the simulated value 
of the particle position: x(nr) = X)j=i(=t)j Pj/2- 
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The sample paths of x{nr) are the piecewise constant random functions of time 
which are defined on the intervals \jr, jr + r) (0 < j < n) and satisfy the condition 
x(0) = (see figure 1(a)). The behavior of these sample paths at r — > is determined 
by two major factors. The first consists in the absence of finite fractional moments of 
pj. This means that the increments Ar^jr) have no characteristic scale and arbitrary 
large values of Ar/(jr) are in principle possible. The second is that the scale function 
a(r) tends to zero very rapidly as r decreases. Since A^(jr) oc a(r), this implies that 
for small enough r the difference between the increments in neighboring intervals is in 
general negligible. Therefore, one may expect (and this is confirmed by our simulations) 
that in the limit r — > the competition of these factors leads to the condition x(t) = 
with t < t tr , where t tT is a random time in which the absolute value of the increment 
becomes bigger than any preassigned positive number (see figure 1(b)). We interpret 
the discontinuous jump of x(t) which occurs at t = ttr as the transition of a particle to 
the absorbing state. 
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Figure 1. Illustrative examples of the sample paths of the particle position x(t) for a 
finite time increment r (a) and in the limit r —¥ (b). The vertical arrow in figure 1(a) 
means that the jump of x(t) at t = At exceeds the figure frame. In contrast, the arrows 
in figure 1(b) indicate that at t — t tr an infinite jump of x(t) occurs in the positive 
or negative direction of the axis x. The horizontal line with x(t) = and t < t tI 
corresponds to the surviving state, and the vertical line to the absorbing state. 



Illustrative histograms of the probability distribution of particle positions are shown 
in figure 2 for two different values of the time increment r. The height of the mth bin 
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Figure 2. Histograms of the probability distribution of particle positions obtained by 
numerical simulations of the Langevin equation (TT|). The parameters of the simulation 
are t = 1, r = 1, s = 2, N = 10 6 , I = 5, e = 0.2, and r = 0.25 (a) and t = 10~ 3 (b). 
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Figure 3. Time dependence of the survival probability. The solid line represents the 
theoretical dependence s -4 / 1 " (for r = 1, s = 2), and the circular symbols indicate the 
results of the numerical simulation. The latter are obtained by calculating the height 
of the central bin in figure 1(b) for different values of the dimensionless time t. 

(m = 0, ±1, . . . , ±(Z — 1)) represents the probability W m to find a particle in this bin 
at time t = nr. The probability W m is determined as the ratio N m /N, where N is the 
overall number of numerical experiments, i.e. calculations of x(t), and N m is the number 
of those experiments in which x(t) falls into this bin, i.e. when x(t) G (me— e/2, me+e/2) 
with e being the bin width. In contrast, the heights of the leftmost and rightmost bins 
represent the probabilities W-i and W\ that x(t) < —le + e/2 and x(t) > le — e/2, 
respectively. According to these definitions Y} m =-i W m = 1- The probabilities W m 
(m 7^ 0, ±Z) rapidly tend to zero as r decreases, and for small enough r the histogram 
takes the form with only three visible bins (e.g., in figure 2(b) max W m ~ 10~ 4 , Wq ~ 0.5, 
W-i « Wi w 0.25). Hence, the less is r the better the condition Wo + W-i + Wi — 1 
holds. The probability Wo approaches the survival probability s~'/ r (see figure 3), and 
its independence on e (if e is not too small for a fixed r) yields W = Pl e dxV(x, t), i.e. 
the theoretical result ( 1231) indeed holds. Finally, the facts that W m (m ^ 0, ±Z) approach 
zero as r decreases and W_i + Wi ~ 1 — s~ t l r (if I is not too large for a fixed r) capture 
the main properties of A(x,t), namely, A(x,t) — > as k — > and Ao(t) = 1 — s~ t ' r . 
Thus, the simulation results are in full agreement with our theoretical predictions. 

5. Conclusion 

In summary, we have incorporated a class of super-heavy-tailed noises into the Langevin 
equation formalism. These noises arise from the super-heavy-tailed probability densities 
whose all fractional moments are infinite. Due to this feature, super-heavy-tailed noises 
induce two probabilistic states of a particle, survived and absorbed ones. Thus the 
Langevin equation driven by these noises becomes a useful tool for studying a number of 
physical phenomena such as absorption, evaporation, annihilation, etc. We have solved 
analytically the corresponding generalized Fokker-Planck equation and have confirmed 
our theoretical predictions via the numerical simulation of this Langevin equation. 
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